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^ \ Abstract 

In this pair of papers (Part I and Part II in this issue), we investigate the issue of power control 
and subcarrier assignment in a sectorized two-cell downlink OFDMA system impaired by multicell 
interference. As recommended for WiMAX, we assume that the first part of the available bandwidth is 
likely to be reused by different base stations (and is thus subject to multicell interference) and that the 
second part of the bandwidth is shared in an orthogonal way between the different base stations (and is 
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5/5 ' thus protected from multicell interference). 
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Although the problem of multicell resource allocation is nonconvex in this scenario, we provide in 

£f) . Part I the general form of the global solution. In particular, the optimal resource allocation turns out to 
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qq ■ be "binary" in the sense that, except for at most one pivot-user in each cell, any user receives data either 
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in the reused bandwidth or in the protected bandwidth, but not in both. The determination of the optimal 
resource allocation essentially reduces to the determination of the latter pivot-position. 
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OFDMA Networks, Multicell Resource Allocation, Distributed Resource Allocation. 

I. Introduction 

We consider the problem of resource allocation in the downlink of a sectorized two-cell OFDMA 
system with incomplete Channel State Information (CSI) at the Base Station (BS) side. In principle, 
performing resource allocation for cellular OFDMA systems requires to solve the problem of power and 
subcarrier allocation jointly in all the considered cells, taking into consideration the interaction between 
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users of different cells via the multicell interference. Unfortunately, in most of the practical cases, this 
global optimization problem is not convex and does not have, therefore, simple closed-form solution. 
Practical alternative methods must thus to be proposed to perform the resource allocation. Most of the 
works in the literature on multicell resource allocation assumed perfect CSI on the transmitters side. 
In flat-fading scenarios with multi-user interference, a number of interesting alternative methods have 
been proposed in the literature. One of them is the geometric programming (GP) approach proposed 
in [1] for centralized power control scenarios. The author of this work showed that at high SNR, the 
GP technique turns the nonconvex constrained optimization problem of power control into a convex, 
thus tractable, optimization problem. Another efficient resource allocation technique was proposed in 
[2] for decentralized power control scenarios. This technique is based on a min-max formulation of 
the optimization problem, and is adapted to ad-hoc networks contexts. Unfortunately, the two above 
mentioned techniques are mainly intended for flat-fading scenarios, and are not directly suitable to general 
cellular OFDMA contexts. To the best of our knowledge, only few works investigate OFDMA multicell 
resource allocation. Authors of [3] addressed the optimization of the sum rate performance in a multicell 
network in order to perform power control and user scheduling. In this context, the authors proposed 
a decentralized algorithm that maximizes an upperbound on the network sum rate. Interestingly, this 
upperbound is proved to be tight in the asymptotic regime when the number of users per cell is allowed 
to grow to infinity. However, the proposed algorithm does not guaranty fairness among the different users. 
In [4], a centralized iterative allocation scheme allowing to adjust the the number of cells reusing each 
subcarrier was presented. The proposed algorithm does not suppose the so called "reuse partitioning" 
scheme but nonetheless it promotes allocating subcarriers with low reuse factors to users with bad channel 
conditions. It also provides an interference limitation procedure in order to reduce the number of users 
whose rate requirements is unsatisfied. Authors of [5] considered the problem of subcarrier assignment 
and power control that minimize the percentage of unsatisfied users under rate and power constraints. 
For that sake, a centralized algorithm based on reuse partitioning was proposed. In this algorithm, the 
reuse factor of the far users next to the cell borders is adapted according to the QoS requirements and the 
problem parameters. Other dynamic resource allocation schemes were proposed in [6]-[10]. The authors 
of [9] and [10] have particularly discussed the issue of frequency reuse planning. It is worth mentioning 
here that neither of the above cited works [4]-[10] provided analytical study of the performance of their 
respective proposed schemes. The issue of power control in distributed cooperative OFDMA networks 
was addressed in [11]. However, the proposed solution assumes that subcarrier allocation is performed 
independently from the power control. The solution is thus suboptimal for the problem of resource 
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allocation for OFDM A networks, and a general solution for both power control and frequency resource 
allocation remains to be provided. 

In contrast to previous works where perfect CSI was assumed, authors of [12] assumed the knowledge 
of only the statistics of users' channels and proposed an iterative algorithm for resource allocation in the 
multicell context. In this algorithm a frequency (or subcarrier) reuse factor equal to one was chosen, which 
means that each cell is supposed to use all available subcarriers. This assumption relatively simplifies 
solving the problem of multicell OFDMA resource allocation. A similar iterative multicell allocation 
algorithm was proposed in [13] and its convergence to the optimal solution of the multicell resource 
allocation problem was proved based on the framework developed in [14]. 

In this paper, our aim is to characterize the resource allocation strategy (power control and subcarrier 
assignment scheme) allowing to satisfy all users' rate requirements while spending the least power at the 
transmitters' side. Similarly to [12], we investigate the case where the transmitter CSI is limited to some 
channel statistics. However, contrary to [12] which assumes a frequency reuse factor equal to one, our 
model assumes that a certain part of the available bandwidth is shared orthogonally between the adjacent 
base stations (and is thus "protected" from multicell interference) while the remaining part is reused by 
different base stations (and is thus subject to multicell interference). Note that this so-called fractional 
frequency reuse is recommended in a number of standards e.g. in [15] for IEEE 802.16 (WiMax) [16]. A 
similar reuse scheme is adopted in the recent work [17] which addresses the problem of power allocation 
in a 2-cell OFDMA system in order to maximize the system sum rate under a total power constraint. The 
method proposed by the authors of [17] to tackle the latter problem is based on a game theory approach 
and it assumes that subcarrier assignment is fixed in advance. 

As opposed to [17], our work considers the problem of joint optimization of power allocation and 
subcarrier assignment under the aforementioned frequency reuse scheme. We also assume that each user 
is likely to modulate in each of the two parts of the bandwidth (the protected and the non protected 
parts). Thus, we stress the fact that i) no user is forced to modulate in a single frequency band, ii) we 
do not assume a priori a geographical separation of users modulating in the two different bands. On the 
opposite, we shall demonstrate that such a geographical separation is actually optimal w.r.t. our resource 
allocation problem. In this context, we provide an algorithm that permits to compute the optimal resource 
allocation. 

The paper is organized as follows. In Section II we present the system model. In Section III we 
consider the problem of resource allocation in a single cell assuming that the interference generated by 
the other cells of the network is fixed. The problem consists in minimizing the transmit power of the 
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considered cell assuming a fixed level of interference such that the rate requirements of users of this cell 
are satisfied and such that the interference produced by the cell itself is less than a certain value. Although 
resource allocation for users of the network requires in general solving a multicell optimization problem, 
the single cell problem of Section III turns out to be a useful tool to solve the more complicated multicell 
problem. Theorem 1 gives the solution to this single cell optimization problem. Except for at most one 
"pivot" user in the considered cell, any user receives data either in the interference bandwidth or in the 
protected bandwidth, but not in both. In Section IV we introduce the joint multicell resource allocation 
problem. This problem is equivalent to jointly determining the resource allocation parameters of users 
belonging to different interfering cells, such that all users' rate requirements are satisfied and such that 
the total transmit power is minimized. Theorem 2 characterizes the solution to this optimization problem 
as function of a small number of unknown parameters. The solution turns out to have in each cell the 
same binary form as the solution to the single cell problem. Although this geographical separation is 
frequently used in practice, no existing works prove the optimality of such a scheme to our knowledge. 
Subsection IV-C provides a method to calculate the optimal resource allocation. Finally, Section V is 
devoted to the numerical results. 

II. System Model 

A. OFDMA Signal Model 

We consider a downlink OFDMA sectorized cellular network. In order to simplify the presentation of 
our results, the network is supposed to be one-dimensional (linear) as in a number of existing studies 
[12], [18], [19], [20], [21]. The motivation behind our choice of the one-dimensional network is that such 
a simple model can provide a good understanding on the problem while still grasping the main aspects of 
a real-world cellular system. It provides also some interesting guidelines that help to implement practical 
cellular systems. Generalization to 2D-networks is however possible (though much more involved) and 
is addressed in a separate work [22]. We consider the case of sectorized networks i.e., users belonging to 
different sectors of the same cell are spatially orthogonal [23]. In this case, it is reasonable to assume that 
a given user is only subject to interference from the nearest interfering base station. Thus, we focus on 
two interfering sectors of two adjacent cells, say Cell A and Cell B, as illustrated by Figure 1. Denote by 
D the radius of each cell which is assumed to be identical for all cells without restriction. We denote by 
K A and K B the number of users in Cell A and B respectively. We denote by K = K A + K B the total 
number of users in both cells. Each base station provides information to all its users following a OFDMA 
scheme. The total number of available subcarriers is denoted by N. For a given user k € 1,2,..., K A 
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Figure 1. Two-Cell System model 

in Cell A, we denote by the set of indices corresponding to the subcaniers modulated by k. yS k is a 
subset of {0, 1, . . . , N — 1}. By definition of OFDMA, two distinct users k, k' belonging to Cell A are 
such that Nfc Pi = 0. For each user k € {1, . . . , K A } of Cell A, the signal received by k at the nth 
subcarrier (n G DM^) and at the mth OFDM block is given by 



where s k (n, to) represents the data symbol transmitted by Base Station A. Process w k (n, m) is an additive 
noise which encompasses the thermal noise and the possible multicell interference. Coefficient H k (n, to) 
is the frequency response of the channel at the subcarrier n and the OFDM block to. Random variables 
H k (n,m) are assumed to be Rayleigh distributed with variance 



Note that the mean value p k does not depend on the subcarrier index. This is satisfied for instance in 
the case of decorrelated channel taps in the time domain. For a given user k, Hk(n,m) are identically 
distributed w.r.t. n, m, but are not supposed to be independent. Channel coefficients are supposed to be 
perfectly known at the receiver side, and unknown at the base station side. However, variances pk are 
supposed to be known at the base station. This type of incomplete CSI is particularly adapted to fast 
fading scenarios. In such a context, sending feedback containing the instantaneous channel gain from 
users to the base station will result in a significant overhead. 



y k {n,m) = H k (n,m)s k {n,m) + w k {n,m), 



(1) 



Pk = E[\H k (m,n)\ 2 }. 



(2) 
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As usual, we assume that p k vanishes with the distance between Base Station A and user k, based 
on a given path loss model. In the sequel, it is convenient to assume (without restriction) that users 
k = 1,2, . . . ,K A are numbered from the nearest to the base station to the farthest. Therefore, for all 
users k in Cell A, 

Pi> P2> ■■■> Pk a ■ (3) 

B. Frequency Reuse 

The frequency reuse scheme is illustrated by Figure 1. In practical cellular OFDMA systems, it is 
usually assumed that certain subcarriers n G {0, . . . N — 1} used by Base Station A are reused by the 
adjacent Cell B. Denote by 3 this set of "Jnterfering" subcarriers, 3 C {0, . . . , N — 1}. If user k modulates 
such a subcarrier n G 3, the additive noise Wk(n,m) contains both thermal noise of variance a 2 and 
interference. Therefore, the variance of Wk(n,m) depends on k and is crucially related to the position 
of user k. We thus define for all n G 3 

E[\w k (n,k)\ 2 ] =a\. 

Note that a\ is assumed to be a constant w.r.t. the subcarrier index n. This assumption is valid in 
OFDMA multicell systems using frequency hopping or random subcarrier assignment as in WiMax. If 
users k = 1,2... K A are numbered from the nearest to the base station to the farthest, it is reasonable 
to assume that 

a\ < a\ < . . . < a\ A , (4) 

meaning that the farthest users experience more multicell interference. The reuse factor a is defined as 
the ratio between the number of reused subcarriers and the total number of available subcarriers: 

card(J) 

a = ^T 

so that 3 contains aN subcarriers. The remaining (1 — a)N subcarriers are shared by the two cells, A 
and B, in an orthogonal way. We assume that of these subcarriers are used by Base Station A 

only and are forbidden for B. Denote by T ^ this set of "^Protected" subcarriers. If user k modulates 
such a subcarrier n G 7 a, the additive noise Wk(n,m) contains only thermal noise. In other words, 
subcarrier n does not suffer from multicell interference. Then we simply write E[|^(n, m)\ 2 } = a 2 , 
where a 2 is the variance of the thermal noise only. Similarly, we denote by Tb the remaining 
subcarriers, such that each subcarrier n G T_b is only used by Base Station B, and is not used by A. 
Finally, 3 U Ta U 3^b = {0, . . . , N — 1}. Moreover, let g^ i (resp. 5^2) t> e the channel Gain to Noise 
Ratio (GNR) in band 3 (resp. 9a), namely # M = p k /a 2 k (resp. g K2 = Pk/o~ 2 ). 
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C. Resource Allocation Parameters 

Of course, for a given user k of Cell A, the noise variance cr? depends on the particular resource 
allocation used in the adjacent Cell B. We assume that a\ is known at Base Station A, and that a given 
user may use subcarriers in both the "interference" bandwidth 3 and the "protected" bandwidth T a- We 
denote by j^N (resp. 7^f 2 A0 tne number of subcarriers modulated by user k in the set 3 (resp. 7a)- 
In other words, 

Tjji = card(J n X fe )/iV 7 £ 2 = card^ n N fc )/JV. 

Note that by definition of 7^ and 7^5 2 , Sfc7£fi ^ a an ^ J2k1k2 — an ^ tnat tne superscript 

^4 (or i?) is used to designate the cell in which user k is located. We assume in the sequel without 
restriction that the sharing factors {7^1,7^2}* WQ continuous real-valued variables and can take on any 
value in the interval [0, 1]. Furthermore, we assume that a given user k of Cell A can modulate in both 
bands 3 and 7a using distinct powers in each band. For any modulated subcarrier n G N^, we define 
Pfa = E[\s k (n, m)\ 2 } if n £ 3, P^ 2 = E[\s k (n, m)\ 2 } if n € 7 a- Similarly, denote by = ^,i P k,i the 
average power transmitted to user k in 3 if % = 1 and in 7a if i = 2. "Setting a resource allocation for Cell 
A" means setting a value for parameters {7^, 7^ ^kii Pk2\k=i...K A , or equivalently for parameters 

bk,i^k,2^ v w^ 2 } k=1 ... K A. 

Remark 1. As we stated above, the sharing factors 7^, 7^ 2 are assumed in our model to be real 
numbers. This assumption does not necessarily contradict the fact that each user can be assigned only 
integer number of subcarriers during the transmission of each OFDM symbol. Indeed, once the real-valued 
{7^, 7^ 2 }fc are determined, the practical subcarrier assignment can be done in several ways [12]. One 
possible way consists in allocating subcarriers to users according to some frequency hopping pattern. 
In this case, the specific subset of subcarriers assigned to each user varies from one OFDM symbol to 
another in such a way that the average number of subcarriers modulated by each user k in bands 3 and 
7a is equal to and 7^ 2 ^ respectively. The latter frequency-hopping-based subcarrier assignment 

scheme is assumed in this paper. 

D. Multicell Interference Model 

We define now more clearly the way interference levels a 2 ,... , o 2 k a depend on the adjacent Base 
Station B. In OFDMA system models which assume frequency hopping like Flash-OFDM system ([23] 
Chapter 4, page 179-180, [24]), it is straightforward to show that for a given user k of Cell A, interference 
power a 2 does not depend on the particular resource allocation in Cell B but only on i) the position of 
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user k and ii) the average power Qf = Yl?k=i ^ki transmitted by Base Station B in the interference 
bandwidth J. More precisely, 

f \H k (n,m)\ 2 }Qf + a 2 (5) 



a 2 k =E 



where Hk(n,m) represents the channel between Base Station B and user k of Cell A at frequency n 
and OFDM block m. In particular, E \Hk(n,m)\ 2 only depends on the position of user k and on the 
path-loss exponent. 

III. Single Cell Resource Allocation 

Before tackling the problem of joint optimal resource allocation in the two considered cells, it is useful 
to consider first the simpler single cell problem. The single cell formulation focuses on resource allocation 
in one cell, and assumes that the resource allocation parameters of users in the other cell are fixed. 

A. Single Cell Optimization Problem 

Assume that each user k has a rate requirement of Rk nats/s/Hz. Our aim is to optimize the resource 
allocation for Cell A which i) allows to satisfy all target rates Rk of all users, and ii) minimizes the 
power used by Base Station A in order to achieve these rates. Considering a fast fading context (i.e. 
channel coefficients Hk(n, m) vary w.r.t. m all along the code word), we assume as usual that successful 
transmission at rate Rk is possible provided that Rk < Ck, where Ck denotes the ergodic capacity 
associated with user k. Unfortunately, the exact expression of the ergodic capacity is difficult to obtain 
in our context due to the fact that the noise-plus-interference process (wk(n,m)) n ^ m is not a Gaussian 
process in general. Nonetheless, if we endow the input symbols Sk(n,m) with Gaussian distribution, the 
mutual information between Sk(n, m) and the received signal yk{n, m) in equation (1) is minimum when 
the interference-plus-noise Wk(n, m) is Gaussian distributed. Therefore, the approximation of the multicell 
interference as a Gaussian random variable is widely used in the literature on OFDMA (see for instance 
[12], [25], [26]) as it provides a lower bound on the mutual information. For a given user k in Cell A, the 
ergodic capacity in the whole band is equal to the sum of the ergodic capacities corresponding to both 
bands J and "Pa- For instance, the part of the capacity corresponding to the protected band 7 'a is equal to 
7^ 2 E log (l + -P^ 2 |g *y 2 ) =7^2 E log(l+ W ^ |gfc( ^' m)|2 ) , where factor 7^ 2 traduces the fact 
that the capacity increases with the number of subcarriers which are modulated by user k. In the latter 
expression, the expectation is calculated with respect to random variable \ H >°(™> n )\ _ Now, gfc (™' n ) h as 
the same distribution as = gk^Z, where Z is a standard exponentially distributed random variable. 
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log (l + g k ,J^Z 



K-" 



+ 7j£ 2 E 



w k A 2 



log \l+g K2 -^-Z 



(6) 



Finally, the ergodic capacity in the whole bandwidth is equal to 

Cfc(7M, lk,2, <i, < 2 ) = 7m E 

The quantity Q A defined by 

Q A = £«i + <2) (7) 
fc=i 

denotes the average power spent by Base Station A during one OFDM block. The optimal resource 
allocation problem for Cell A consists in characterizing {7^, 7^ 2 , W k A 2 \ k =\...K A allowing to satisfy 
all rate requirements of all users {R k < C k ) so that the power Q A to be spent is minimum. Furthermore, 
as we are targeting a multicell interference scenario, it is also legitimate to limit the interference which 
is produced by Base Station A. Therefore, we introduce the following "low nuisance constraint": The 
power Q A = J2k ^ki which is transmitted by Base Station A in the interference band 3 should not 
exceed a certain nuisance level Q, which is assumed to be a predefined constant imposed by the system's 
requirements. The introduction of this constraint will be later revealed useful in Section IV when studying 
the solution to the joint multicell resource allocation problem. The single cell optimization problem can 
be formulated as follows. 

Problem 1. Minimize Q A w.r.t. {l k i,l k 2->W kl ,W A 2 }k=i...K A under the following constraints. 
CI : Vk,R k <C k C4 : 7^ > 0, 7 £ 2 > 

K A 

C2 : = « C5 : > 0, W A 2 > 0. 

k=i 

C3:^ 7 ^ = i T ^ C6:f>&<0. 

k=l k=l 

Here, CI is the rate constraint, C2-C3 are the bandwidth constraints, C4-C5 are the positivity 
constraints. Note that C6 is the low nuisance constraint imposed only on the power transmitted in the 
non protected band 3. The particular case where the maximum admissible nuisance level is set to Q = +00 
would correspond to a "selfish" resource allocation: Base Station A may transmit as much power as needed 
in the interference band 3 without caring about the nuisance which it produces on the adjacent cell. Note 
that in Problem 1 no power constraint is imposed on the total power Q A transmitted by the base station in 
the two bands. Note also that the constraint set (the set of all feasible points) associated with Problem 1 is 
not empty as it contains at least the following trivial solution. This trivial solution consists in assigning zero 
power W A X = on the subcarriers of the non protected band J (so that constraint C6 will be satisfied), 
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and in performing resource allocation only in the protected band "Pa- The main reason for expressing the 
resource allocation problems in terms of parameters l^i^ki (i = 1,2) instead of ^v^ki * s tnat tne 
ergodic capacity = ^(7^, W^, 7j^ 2 > ^^2) * s a concave function of 7/fj,W / fcV As a consequence, 
the constraint set is convex and Problem 1 is a convex optimization problem in {7^, 7^ 2 , Wj^ 2 }k- 
Obviously, finding the optimal parameter set {7^, 7^ 2 , W£ 2 }k is equivalent to finding the optimal 
{^i, 7j£ 2 . P k,v p k,2}k thanks to the simple relation = J^P^, i = 1,2. 

B. Optimal Single Cell Resource Allocation 

In order to solve convex Problem 1, we use the Lagrange Karush-Kuhn-Tucker (KKT) conditions. 
Define the following function on R + 



E 



1+xZ 



(8) 



It can be shown that function /(x) is increasing from to 00 on E + . The following theorem provides 
the general form of any global solution to Problem 1. Its proof is provided in Appendix A. 

Theorem 1. Any global solution {7^, 7^ 2 , ^kv W k A 2 \k=i...K A to Problem 1 is as follows. There exists 
an integer L G {1, . . . , K A } and three nonnegative numbers 0i, /3 2 and £ such that 
1) For each k < L, 



_ -1 f -l / 9k,l R 

k,i ~ 9 k J ( 



2) For eac/i k > L, 



7fc,i 



1+f 



E 



76 = 



\og{\ + g K1 P^Z 

Pk,2 = 9klr l (9k^2) 

Rk 



Pk A 2 = 



(9) 



lk,2 



(10) 



E 



log 1 + gt 



3) For k = L 



Pi 



-1 e-l I 9k,l ,, 
k,i = 9 k J ( TTTfo. 



k-i 



1+r 



r,2 = 9k~ir 1 (9k,2P2) 

K A 



ll2 



1 - a 



(11) 



z=fc+i 



where 0\, /3 2 a«<i £ are the Lagrange multipliers associated with constraints C2, C3 and C6 respectively. 
Determination of L, (5\, /3 2 and £ is provided by Proposition 1. 



Comments on Theorem 1: 
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a) Theorem 1 states that the optimal resource allocation scheme is "binary": Except for at most one 
user (k = L), any user receives data either in the interference bandwidth J or in the protected 
bandwidth 7 'a, but not in both. Intuitively, it seems clear that users who are the farthest from the 
base station should mainly receive data in the protected bandwidth Ta, as they are subject to an 
significant multicell interference and hence need to be protected. Now, a closer look at our result 
shows that the farthest users should only receive in the protected bandwidth 7a- On the other hand, 
nearest users should only receive in the interference bandwidth 3. 

b) Nonzero resource allocation parameters "l^ii^ki ^ or ^ — ^) and j£ 2 , P£ 2 (for k > L) are 
expressed as functions of three parameters /3i,/?2,£. It can be easily seen from Appendix A that 
/?i,/?2,C are the Lagrange multipliers associated with constraints C2, C3 and C6 respectively. It 
is quite intuitive that, when the admissible nuisance level is large (take for instance Q = +oo), 
constraint C6 holds with strict inequality. Thus, £ = from complementary slackness condition. 
In the general case, the values of parameters f3i, 02, £ can be obtained from KKT conditions. The 
determination of /3i , /?2 , C an d the pivot-user L is addressed in Subsection III-C. 

c) As expected, the optimal resource allocation depends on the resource allocation in Cell B via 
parameters af,... ,ct 2 K a- Joint optimization of the resource allocation in both cells, A and B, is 
investigated in Section IV. 



C. Determination of L, f3\, 02 and £ 
Step 1: Determination of L, fli, [3^ for a fixed value of £. 

To simplify, first assume that the value of Lagrange multiplier £ is fixed. We determine L, 0\, 02 as 
functions of £. Recall that user L is defined as the only user who is likely to modulate in both bands 
J and 7 a- Parameters 7^1,7^2 respectively provide the part of the band J and "Pa which is modulated 
by user L. A first equation is obtained by writing that Cl = Rl i-e-, the rate constraint CI holds with 
equality. Recall that C L is defined by (6) as 7 f ;1 E [log(l + g^P^Z)^ + 7 £ 2 E [log(l + g^Pfe) ■ 
Define for each x > 

C(z)=E[log(l + r 1 (x)Z)]. (12) 



Plugging the expression (1 1) of parameters 7^ ± , 7^ 2 , ^L2 i nto ^ s expression, equality Cl/Rl = 1 
becomes 



a 



E 



Ri 



r< ( 9k 

k<L ^ ^T+£ 



Rl 



+ 



1 - a 



E 

k>L 



Ri 



C{g L ,2fh) 
Rl 



(13) 
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In equation (13), both terms enclosed inside the brackets coincide with 7^ and 7^ 2 respectively. As 
function C(x) is increasing from to 00 on R+, constraints > and 7^ 2 > hold only if 
/?i/(l + £) > a L-i anc ^ 02 > where for each I, af and 6^ the unique positive numbers such that: 

l K A 

\ - Rk _ , v-> Rk 1 — a 

^C( fcl af)- a m ^Oto^)- 2 ■ 

with = = by convention. Note that af is an increasing sequence while bf is a decreasing 
sequence. Furthermore, in order that (13) holds, both (nonnegative) terms should be less than one. Thus, 
a — — t— ^ — r- < and - — — — < 0. As a consequence, f3±/(l + £) < and 

#2 < 6f_i- Finally, 

(iT^ 2 ) G ^-i' a ^ x ^' 6 ^- (14) 

Consider the case where 7^, 7^ 2 are both nonzero, and define the following function for each x > 0: 



It can easily be seen from the KKT conditions derived in Appendix A that 



(15) 



UL,i zF ( IhLfa \ = g Lt2 F(g L M , (16) 



Now using (14) in the above equation along with the fact that F(.) is a decreasing function, one can 
easily see that L can be defined as 

L = min j/ = 1 . . . K A /^F (g^af) < %2 F ( 9l , 2 bf) j . (17) 

In practice, the search for L can be achieved by dichotomy, computing af and bf only for a limited 
number of values of /. Once L is fixed, it is straightforward to show that the system formed by equation 
(16) and (13) admits a unique solution (/?i,/?2)- This is due to the fact that functions C(.) and F(.) are 
monotone. Lagrange multiplier f}\, (5 2 can thus be obtained using classical root search tools. As a remark, 
we note the existence of a rather pathological case, which we do not address in details because of its 
limited importance. To obtain equation (16) we assumed that jf 1 and 7^ 2 are strictly positive. If this is 
not the case, say jf 1 = 0, it turns out that the system (13)-(16) has no solution. However, L can still be 
obtained by (17) and (3\,/3 2 can be easily obtained from (13) which lead to [3\ = (1 + £)°l> P2 = 
For the sake of simplicity, we will still refer to ((3i,/3 2 ) as the unique solution to system (13)-(16), with 
slight language abuse, keeping in mind that we just put /3± = (1 + a L> P2 = ^l-i m tne pathological 
case where such a solution does not exist. This convention will be used throughout the paper without 
restriction. 
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Step 2: Determination of £. 

So far, we proved that for a fixed value of £, the optimal resource allocation is unique and follows 
equations (9), (10) and (11), where L = L(£) is given by (17) and (0i,0%) = (01 (£)>/%(£)) * s tne ur >iq ue 
solution to system (13)-(16). The remaining task is now to determine £. Before addressing this point, it is 
worth providing some insights on the impact of £ or equivalently, on the role of the low nuisance constraint 
C6 on the resource allocation. Recall that £ is the Lagrange multiplier associated with constraint C6. 
From an intuitive point of view, a large value of £ means in some sense that constraint C6 is severely 
restraining, whereas £ = means that constraint C6 has no role and could have been deleted without 
modifying the solution to Problem 1. It turns out that increasing £ has the effect of decreasing the total 
power Qf = ^ki^ki which is transmitted in the interference band. This statement can be proved 
as follows. First, we observe from equation (17) that parameter L = is a non increasing function 
of £. Second, it is straightforward to show that for each k, P kl is a decreasing function of £. Indeed, 
equation (9) implies that it is the composition of an increasing function f~ l (x) and a decreasing function 
£ l— * A (0/(1 +0 (decreasingness of /3i(£)/(l + £) is obtained after some algebra from (13) and (14)). 
Third, W kl = P^Pfc/E log(l + gk^P^Z) is an increasing function of P^V It is thus a decreasing 
function of £ as a composition of an increasing and a decreasing function P£ v Therefore, the presence 
of an active constraint C6 has a double impact on the resource allocation: i) it decreases the number L 
of users who modulate in the interference band J, and ii) it decreases the power Wk,i of each user in 
this band. 

We now determine £. First we propose to compute the resource allocation assuming £ = 0. If the 
corresponding value of Qf is such that Qf < Q, then the procedure stops: KKT conditions are met. 
Otherwise, this means that constraint C6 should be active: £ > 0. From complementary slackness 
condition, C6 should be met with equality : one should determine £ such that = ^k^ki^ki 
coincides with Q: 

k<L 

where Iki^ki WQ defined by (9) and where L = L(£), (3\ = /3i(£)> 02 = 02(0 have been defined 
previously. As mentioned above, Q± is a decreasing function of £ so that the solution £ to equation 
Qi = Q is unique. 

Finally, we conclude that the following proposition holds true. 

Proposition 1. The global solution to the single cell Problem 1 is unique and is given by equations (9)- 
(10)-(11), where parameters L, (3±, fii and £ are unique and determined as follows. 
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1) Assuming £ = 0, evaluate L by (17) and (/3i , /?2 ) as the unique solution to the system of equa- 
tion (13)-(16) satisfying (^,#2) G [a£_ l5 a£] x [6f ,6f_j]. rfcen evaZwate = Ek1k,i P k,v 

2) 5to/7 j/Qi 1 < Q (constraint C6 « mef) otherwise continue. 

3) Evaluate (L, fa, fa, £) as the unique solution to the system of equations (13)-(16)-(17)-(18). 

The above proposition proves that the global solution to the single cell allocation problem is unique 
and provides a procedure to compute it. Algorithm 1 gives a more detailed description of the latter 
procedure and proposes a method to solve the system of equation (13)-(16)-(17)-(18). 

Algorithm 1 Determination of L, fa, fa, £ 
repeat 

L <_ min {l/f± F ( gijia f) < gi ^ F ( gi2b f) } 

(Pi, fa) <- unique solution in [af.^af] x [fef , 6f_ x ] to (13)-(16) 

Qt - SfcJ ^y*,!/- (ftjA) + (« - E£3 s^y) si*/"' (f*ft) 

if Qf > Q then 

Increment £ 
end if 
until Q^ 1 < Q 

return L, fa, fa, ^ 



One still needs to define the way £ should be incremented at the end of each iteration of Algorithm 1 . 
In practice, updating the value of £ can be done by resorting to Newton-like methods which are widely 
used to solve nonlinear equations. 

IV. Joint Multicell Resource Allocation 

A. Optimization Problem 

Our aim now is to jointly optimize the resource allocation for the two cells which i) allows to satisfy 
all target rates Rk of all users, and ii) minimizes the power used by the two base stations in order to 
achieve these rates. The ergodic capacity associated with user k in Cell A is given by equation (6), where 
coefficient 5^1 in that equation coincides with 

9k,l(Ql 



E 



H k (n,m)\ 2 Qf + a 2 
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where H k (n, m) represents the channel between Base Station B and user k of Cell A at frequency n and 
OFDM block m. Coefficient gk,i(Q?) represents user k channel gain to interference -plus-noise ratio in 
the interference band J. Here, gk,i(Qf) not only depends on the position of user k in Cell A, but also 
on the power Qf = Y^k=i ^fi transmitted by the adjacent Base Station B in band J. We now solve 
the following multicell resource allocation problem. 

K c 

Problem 2. Minimize the total power spent by both base stations Q = ^^^(Wk.i + ^£,2) w ^ 

c=A,B k=l 

respect to {7^ v *y% 2 , W£ v W£ 2 } c= a,b under the following constraints. 

' k=l..'.K c 

CI :Vk,R k <C k C4: 7 ^ >0, 7 ^ 2 >0 

C2 : 7m = « C5 : > 0, W^ 2 > 0. 

k=i 

fc=l 

It can be easily seen that the above optimization problem is feasible as soon as a < 1. Indeed, a naive 
but nevertheless feasible point can be easily constructed by forcing each user to modulate in the protected 
band only (force 7^ 1 = for each user). Cells thus become orthogonal, and all users rate requirements 
Rj, can be satisfied provided that enough power is transmitted in the protected band. Unfortunately, the 
ergodic capacity C k of user k is not a convex function with respect to the optimization variables. This is 
due to the fact that the gain-to-noise ratio gk,i(Q\ ) is a function of the resource allocation parameters 
of users belonging to the interfering cell. Therefore, Problem 2 is nonconvex, and cannot be solved by 
classical convex optimization methods. Nonetheless, we manage to characterize its solution. In fact, we 
prove that the solution has the same simple binary form of the single cell optimal solution. 



B. Optimal Resource Allocation 

For each cell c <G {^4, B}, denote by c the adjacent cell (A = B and B = A). The following result is 
proved in Appendix B. 

Theorem 2. 

(A) Any global solution to Problem 2 has the following form. For each Cell c, there exists an integer 
L c G {1, . . . , K c }, and there exist four positive numbers f5\, (3%, £ c , Q\ such that 
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1) For each k < L c , 



pc _„ / n c\-lf-l ( 9k,l(Qj) ac 

n,i - 9k,i{Qi) J \ , tc Pi 



7fc,i 



1 + e 



E 



2) For eac/i > LP, 



7 M 



log (1 + g Kl {Q\)P^Z 



P C K2 = 



7^,2 = 



(19) 



7fc,2 



(20) 



E 



3) For k = L c 



pc 



9k,i{Qi) J 1 i i ^ n Pi 



fe-i 



Pk,2 = 9kir\9k,2P c 2 ) 



7fc,i = « ~ E 7u 



7fc,2 



1 — a 



E ^2- 

l=k+l 



(21) 



(B) For eac/i c = A,B, the system § c (Qf, Qf) formed by the following four equations is satisfied. 
L = mm <! I = 1 . . . K / x ^ Cc F [ -L—— ai ) < g lj2 F (g^ 2 b h 



i + e 



I — : — Pi - gL-,2^ \9L^2P2) 



P c 1 \+l c L c, 2 C(g L , 2 p c 2 ) = R L c 



i + e 

E^,i P M = 3i, 



(22) 
(23) 
(24) 

(25) 



where the values ofj% 1 an J in (25) are the functions of (Pf, (3 2 ,£, c ) defined by equation (19). 
(C) Furthermore, for each c = A,B and for any arbitrary values Qi and Qf, the system of equations 
& c (Qi,Qi) admits at most one solution (L c , f3{, 2 , £ c ). 

Comments on Theorem 2: 

a) The joint multicell resource allocation problem required initially the determination of 4K parameters 
(where K is the total number of users). Theorem 2 allows to reduce the search to only two 
parameters, namely Qf and Qf. Once the value of these parameters is fixed, the resource allocation 
parameters for each user can be obtained from the above results. As a consequence, the only 
remaining task is to determine the value of (Qi,Qf ). This task is addressed in Subsection IV-C. 

b) We observe that the system S c {Q±~,Qf) is very similar to the system obtained in the single cell 
case at equations (13), (16), (17) and (18). In fact, as stated by the proof later, the optimal resource 
allocation in the multicell case can be interpreted as the solution to a certain single-cell problem. 
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c) As a consequence of the above remark, Theorem 2 states that the optimal multicell resource 
allocation scheme has the same "binary" form as in the single cell case. Even if optimal resource 
allocation is achieved jointly for both interfering cells, there still exists a pivot-user L c in each 
Cell c which separates the users modulating respectively in bands J and (P c . 

d) It is worth noticing that this binary resource allocation strategy is already proposed in a number 
of recent standards. One of the contributions introduced by Theorem 2 is the proof that such a 
strategy in not only simple and intuitive, but is also optimal. 

C. Optimal Distributed Algorithm 

Once the relevant values of Q A and Qf have been determined, each base station can easily compute 
the optimal resource allocation based on Theorem 2. As a consequence, the only remaining task is to 
determine the value of (Q A , Qf)- To that end, we propose to perform an exhaustive search on (Q A , Qf)- 

i) For each point (Q A ,Qf) on a certain 2D-grid (whose determination will be discussed later on), each 
base station c = A,B solves the system § C (Q A ,Qf) introduced by Theorem 2. Solving § C (Q A ,Qf) 
for arbitrary values (Q A ,Qf) can be easily achieved by base station c thanks to a simple single-cell 
procedure. Focus for instance on Cell A. 

• Base station A solves the single cell resource allocation Problem 1 assuming that the interference 
level coincides with Qf, and that the nuisance constraint Q is set to Q = Q A . The (unique) solution 
is provided by Theorem 1 and Proposition 1. 

• If the resulting power ^ fc j A iP A i transmitted in the interference band "?a is equal to the nuisance 
constraint Qf (i.e. constraint C6 holds with equality), then the resulting value of (L A , [5 A , /3 A , £ A ) 
coincides with the unique solution to system $ A (Q A , Qf). This claim is the immediate consequence 
of Proposition 1 . 

. If the power Sfc 7jfci-ffc\ * s l ess than Q A (i.e. constraint C6 holds with strict inequality), then 
(L A , f3 A , $2 > £ A ) is clearly not a solution to system B A (Qf,Qf), as equality (25) does not hold. 
In this case, it can easily be seen that $ A (Q A , Qf ) has no solution. The point (Q A ,Qf) cannot 
correspond to a global solution as stated by Theorem 2 and is thus eliminated. 

ii) Base station A evaluates the power 

Q${Qt Qf) = E ^Vm + 

k 

that would be transmitted if the interference level and the nuisance constraint were respectively equal to 
Qf and Q A . This value is then communicated to Base Station B. Base station B proceed in a similar 
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way. 

Hi) The final value of (Qf,Qf) is defined as the argument of the minimum power transmitted by the 
network: 

(Q^,Qf) = arg min Qf(Qf, Qf) + Qf(Qf, Qf) ■ 

(Q?:Q") 

Note that the optimal resource allocation algorithm as described above does not require the intervention 
of a central controlling unit supposed to have access to the two base stations and to users' information 
(position and data rate). We only assume that both base stations can communicate via a special link 
dedicated to this task. The algorithm is thus distributed. This special link will be only used to exchange 
a limited number of messages. Indeed, the only values that need to be exchanged between the two base 
stations are Qf(Qf,Qf) an d Qf(QfiQi) corresponding to the couples (Qf,Qf) for which the two 
systems of equations § A (Qf, Qf) and § B (Qf, Qf) have a solution. Algorithm 2 given below summarizes 
the steps involved in the optimal resource allocation. 

Algorithm 2 Optimal distributed allocation algorithm 



1. Single cell processing 



Cell A 


Cell B 


for each (Qf,Qf) do 


for each (Qf,Qf) do 


bk,n p £ihk - Solve S A (Q?,Qf) using 


bk,v P k,i)i,k <- Solve S B (Qf,Qf) using 


Algorithm 1 with Q = Qf 


Algorithm 1 with Q = Qf 


if Ef=i 7,Vm = Qf then 

Qf(QliQl) <— Ei=l,2 Ejfc=l lk,i P k,i 


ifEf^T^^Qf then 

Qt(Q A i Qf ) ^~ Ei=i,2 Efe=i ik,i p k,i 


end if 


end if 


end for 

2. Cooperation between BS A and B 


end for 


(Qf, Qf) - argmin^B) Qf(Qf, Qf) + Q$(Qf, Qf) 


3. Resource allcoation in each cell 
Cell A 


Cell B 


{lk,vPk,iU - Solve § A (Qf,Qf) using 


{lk,vPk,ihk - Solve § B (Qf,Qf) using 


Algorithm 1 with Q = Qf 


Algorithm 1 with Q = Qf 



Determination of the search domain in (Qf , Qf ). 
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In order to limit the complexity of the proposed approach, the search for (Qi, Qf) should be restricted 
to a certain compact domain, say 

Q\ G [0, Qmax] 

for each c. For instance, a possible value for Q m ax can be defined as the total power that would be spent 
by the two base stations if one would use the naive and suboptimal resource allocation which consists 
in only transmitting in the protected bands 7 'a and Tb- Clearly, the latter value of Q m ax is a constant 
w.r.t. and Qf and can be computed beforehand. A second way to restrict the search domain is to 
make use of a simple suboptimal multicell resource allocation algorithm prior to the use of our algorithm 
(see for instance the suboptimal algorithm defined in Part II of this work). In this case, it is possible to 
restrict the search for (Qi,Qf) to a well-chosen neighborhood of the couple (Qi, Qf ) su bopt provided 
by the suboptimal solution. 
Complexity Analysis. 

In order to get an idea about the cost of applying the optimal allocation, we provide in the following a 
computational complexity analysis of this algorithm as function of the number of users K in the system. 
For that sake, recall that the system of equations $ c (Qf, Qf) must be solved for each possible value of 
(Qi, Qf) inside a 2D-grid contained in a compact interval. 

For each point (Q?,Qf) of the 2D-grid, solving § c (Qf,Qf) can be done by a procedure similar to 
Algorithm 1 . Recall that during each iteration of the latter algorithm, the value of L c should be determined 
by solving L c = min {1 = 1... K° / i+j^F (gi^af) < gipF (giptf) }. This requires that a certain subset 
of parameters {af, bf} should be computed first. It can be shown that the number of operations needed to 
compute af,bf is of order 0(K C ). Furthermore, we argued in Subection III-C that the determination of 
L c can be done by dichotomy, computing a\ and b\ only for a limited number, namely log 2 K c , of values 
of /. The overall complexity of finding L c for a fixed £ c is therefore of the order of 0(K c \og 2 K c ). 

Once L c is determined, the following step of Algorithm 1 consists in solving the system of equa- 
tion (13)-(16) in variables This system of non linear equations can be solved using Newton-like 
iterative methods. One can verify by referring to [27] that the latter system requires a computational 
complexity proportional to 0(K C ). The computational complexity associated with each iteration of 
Algorithm 1 is therefore dominated by the cost of computing L c , which is of order 0(K C log 2 K c ). 
Denote by iVj the number of iterations of Algorithm 1 needed till convergence. We conclude that the 
overall computational complexity of solving §> c (Qi,Qf) is of the order of 0(NiK c log 2 K c ). 

Denote by M the number of couples (Qf, Qf) in the 2D-grid. The overall computational complexity 
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of the algorithm can be obtained by multiplying the cost of solving S C (Q A ,Q B ) for each point of 
the 2D-grid by M the number of points in the grid. The latter overall cost is therefore of the order of 
0{MNiK A log 2 K A )+0{MNiK B log 2 K B ), which is itself of order 0(MNiK log 2 K) in the particular 
case K A ~ K B ~ if/2. 

Note from the above discussion that the determination of the pivot-user L c in each cell for each value 
of (Q A ,Q B ) is one of the costliest operations in solving § c (Q A ,Qf ) and that it dominates the overall 
complexity. This is why we propose in Part II of this work a simplified resource allocation algorithm which 
uses a predefined value for the pivot distance. The simplified algorithm turns out to have a computational 
complexity of order O(K), as opposed to the computational complexity of the optimal algorithm which 
is of the order of 0(MNiK log 2 K). 

V. Simulations 

In our simulations, we considered a Free Space Loss model (FSL) characterized by a path loss exponent 
s = 2 as well as the so-called Okumura-Hata (O-H) model for open areas [28] with a path loss 
exponent s = 3. The carrier frequency is /o = 2AGHz. At this frequency, path loss in dB is given 
by P<ib{x) = 201og 10 (a?) + 100.04 in the case where s = 2, where x is the distance in kilometers 
between the BS and the user. In the case s = 3, pdB{x) = 301og 10 (x) + 97.52. The signal bandwidth B 
is equal to 5 MHz and the thermal noise power spectral density is equal to iVo = —170 dBm/Hz. Each 
cell has a radius D = 500m and contains the same number of randomly distributed users (K A = K B ). 
The rate requirement of user k in bits/sec/Hz is designated by R^. The distance separating each user 
from the base station is considered a random variable with a uniform distribution on the interval [0,D]. 
The joint resource allocation problem for Cells A and B (Problem 2) was solved for a large number 
of realizations of this random distribution of users and the values of the resulting transmit power were 
averaged. Computing the mean value of the total transmit power with respect to the random positions of 
users is intended to get results that do not depend on the particular position of each user in the cell but 
rather on global information about the geographic distribution of users in the cell. We give now more 
details on the way simulation were carried out. 

Define x as the vector containing the positions of all the users in the system i.e, x = (xi,X2,..., 
xkc)c=a,b- Recall that V/c, Xk is a random variable with a uniform distribution on [0, D]. For each real- 
ization of x, denote by Qr(x, a) the minimal total transmit power that results from a global solution to the 
multicell resource allocation problem (Problem 2) i.e., Q T (x, a) = E C =A,B (E£i W k,i + Ek=L^ W k,^) 
where (7^ 1 , W^,^ 2 , W£ 2 ) c( z{^B},k=i, ...,K<= i s a global solution to Problem 2 described by Theorem 2. 
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Define r t = ^2k=i RkB as the sum rate of the users of Cell c measured in bits/sec. We consider first 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

a a 



Figure 2. Power vs. a for s = 2, D = 500 m, K A = Figure 3. Power vs. a for s = 2, D = 500 m, K A = 
K B = 25, rt = 5 Mbps K B = 25, n = 10 Mbps 

the case where all the users have the same rate requirement R\ = R% = ■ ■ ■ = Rk c - Figures 2 and 3 
represent, for a sum rate requirement of r t = 5 Mbps (Mega bits/sec) and r t = 10 Mbps respectively 
and assuming s = 2, the mean value of Qt(x, a) normalized by its minimum value w.r.t a i.e., the ratio 
E x [Qt(x, a)]/E x [Qr(x, «o)]> where «o is the value of the reuse factor a that minimizes E[Qr(x, a)]. 
Figures 4 and 5 plot the same quantity for r t = 5 Mbps and rt = 10 Mbps respectively, but with the 
difference that it assumes s = 3. The error bars in the aforementioned four figures represent the variance 
of Q r (x,a) i.e., E x [(Q T (x, a) - E x [Q T (x, a)]) 2 ]. 

For each value of x and of the reuse factor a, Qt(x, a) was computed using the optimal resource 
allocation algorithm of Section IV. Power gains are considerable compared to the extreme cases a = 
(the available bandwidth is shared in an orthogonal way between Cells A and B) and a = 1 (all the 
available bandwidth is reused in the two cells). Note also that for r t = 10 Mbps, qo the optimal value of 
the reuse factor that minimizes Qr(x, a) is smaller than the optimal value of the reuse factor for r = 5 
Mbps. This result is expected, given that higher values of r t will lead to higher transmit powers in order 
to satisfy users' rate requirements, and consequently to higher levels of interference. More users will 
need thus to be protected from the higher interference. For that purpose, a larger part of the available 
bandwidth must be reserved for the protected bands 7 'a and 7b- We also remark that in the case where 
s = 3, the value of the reuse factor cko is larger than its value for s = 2. This is due to the fact that 
when the path loss exponent is larger, the interference produced by the adjacent base station will undergo 
more fading than in the case when the path loss exponent is smaller. As a result, less users need to be 
protected from interference in the case s = 3 compared to the case s = 2. (see table V which provides, 
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S=3, D=500m, K A =K B =25 s=3, D=500m, K A =K B =25 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

a a 



Figure 4. Power vs. a for s — 3, D = 500 m, K A = Figure 5. Power vs. a for s = 3, D — 500 m, K A — 
K B = 25, r t = 5 Mbps K B = 25, n = 10 Mbps 

in the two cases, the percentage of protected users to the total number of users for r t = 5 and r t = 10 
Mbps, provided that the corresponding value of «o is used in each case). 





s = 2 


s = 3 


r t — 5 Mbps 


19.8% 


11.6% 


rt = 10 Mbps 


30.0% 


18.7% 



Table I 

Percentage of the protected users to the total number of users 



We now compare the performance of our proposed resource allocation with the distributed scheme 
proposed in [13]. The latter scheme assumes a reuse factor a equal to one (all the subcarriers can be 
reused in all the cells), in contrast to our scheme which uses an optimized value of a. Figure 6 plots the 
average total transmit power E[Q^5^(X, oq)] that results when our proposed scheme is applied compared 
to the power that results from applying the scheme of [13]. This comparison was carried out assuming 
K A = K B = 25, s = 2 and r t = 5 Mbps. The gain obtained when the proposed scheme is applied is 
clear from the figure, and it increases with respect to the total rate r t . We consider now the case when 
the rate requirement is not the same for all users. In particular, we assume that the rate requirement 
of each user is a random variable that can take on one of two values with the same probability. For 
example, consider the case K A = K B = 25 and assume that the rate requirement of each user can either 
be equal to 250 kbps (kilo bits/sec) with probability 0.5 or to 150 kbps with the same probability. This 
means that the mean rate per user is equal to 200 kbps and that the mean total rate per sector is equal 
to r*t = 25 * 200 kbps = 5 Mbps. Figure 7 represents, assuming s = 2, the mean value of Qt(x.,o) 
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Figure 6. Comparison between the proposed optimal scheme Figure 7. Power vs. a for s = 2, D — 500 m, K A 
and the distributed scheme of [13] for K A = K B = 25. K B = 25 assuming random rate requirements. 



normalized by its minimum value w.r.t a i.e., the ratio E x [C}t(x, a)]/E x [Qj'(x, «o)]> where «o is the 
value of the reuse factor a that minimizes E[Qj-(x, a)]. The error bars in the above figure represent the 
variance of Qt(x, a) i.e., E x [(Qy(x, a) — E x [C}t(x, a)]) 2 ]- By comparing Figures 2 and 7 we note that 
the normalized mean value E x [Qr(x, a)]/E x [Qr(x, ao)] is practically the same in the two figures. Only 
the variance E x [(Qy(x, a) — ExfQrlX) a)]) 2 ] is slightly different (its value is slightly larger in the case 
of random rate requirements). 

VI. Conclusions 

In this paper, the resource allocation problem for a sectorized downlink OFDMA system has been 
studied in the context of a partial reuse factor a G [0,1]. The general solution to the (nonconvex) 
optimization problem has been provided. It has been proved that the solution admits a simple form and 
that the initial tedious problem reduces to the identification of a restricted number of parameters. As 
a noticeable property, it has been proved that the optimal resource allocation policy is "binary": there 
exists a pivot-distance to the base station such that users who are farther than this distance should only 
modulate protected subcarriers, while closest users should only modulate reused subcarriers. 



Appendix A 
Proof of Theorem 1 



When the resource allocation parameters of users in Cell B are fixed, it is straightforward to show that 
the ergodic capacity Ck = Cfc(7j^x' 7j^2> ^kf' ^k2) defined by (6) is a concave function of j£ v j£ 2 , 
W^i, (and hence — Cj^^ x , 7j(f 2 , Wfcu 2 ) is convex). This is essentially due to the fact that 
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9k,i = 9k,i{Qf) can be treated as a constant and does not depend on the optimization parameters. Thus, 
the single cell resource allocation problem (Problem 1) is convex in {7^, 7^ 2 > Wj^ 2 }ke{i,...,K A }- 
In the following, we derive the KKT conditions in order to obtain the general form of the solution and 
to prove the existence of L, f3\, 02, £ as stated by Theorem 1. In particular, we prove that any optimal 
resource allocation is binary i.e., there exists a certain pivot-integer L such that 7^ 2 = for k < L 
and 7^ = for k > L. Furthermore, we prove that there exist three parameters f3\, 02 and £ such that 
equations (9), (10) and (11) hold. As explained above, /3i,/?2,£ are the Lagrange multipliers associated 
with constraints C2, C3 and C6 respectively. 
KKT Conditions for Problem 1 

In order to simplify the notations and since we are only interested in users of Cell A, we sim- 
ply omit the superscript A in the sequel and define Q = Q A , 7^ = 7^, etc. Denote by the 
vector of resource allocation parameters of users in Cell A i.e, = [(W) T , (7) T ] T where W = 
[Wi,i, W 1>2 , W K A yl , W k a )2 ] t and 7 = [71,1, 71,2, • • • , 1k a ,i,1k a ,2] T ■ The associated Lagrangian is 



equal to: 



K A / 
k=l V k 



^^fe,i7fc,i - ^^,27fc,2 - - ^2vk,2W k ,2 +£^W fc) i. (26) 

Pc Ic Jc Ic \\j 

where (3\, and £ are the Lagrange multipliers associated respectively with constraints CI, C2, C3 
and C6 of Problem 1, and where 1^1, 1^2, /"fc,i> /"fc,2 are the the Lagrange multipliers associated with the 
positivity constraints of 7^1,7^2, Wfc,ij W^p respectively. In the expression of Cf~, a technical difficulty 
arises from the fact that function 7fc j jE log ^1 + g^^^^-Z^j is not differentiable at point 7^ = 0. One 
can easily overcome this issue by replacing the non-negativity constraint 7^ > by the strict positivity 
constraint 7^ > eo, for an arbitrary eo > 0. However, as this point is essentially technical, we simply 
put eo = with slight lack of rigor. This assumption will simplify the presentation without changing the 
results. The complete proof that does not make this simplifying assumption can be found in [29]. We 
now apply the Lagrange-Karush-Kuhn-Tucker conditions to characterize the optimal vector x^. Taking 
the derivative of (26) with respect to Wk,i and 7^ (i = 1, 2) leads to 

Z 



1 - h9k,i^ 



log 1 + g k , 



W fc , 

' lk,i 



z\- 



l + 9k,^fZ 
1 + 9*^2 



- Vk,i + & = 
+ A - f/c,i = 



(27) 
(28) 
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where Si = 1 if i = 1 and <5j = if i = 2. We can easily show that the constraint i?£ < C k must hold with 
equality, and is always active in the sense that the Lagrange multiplier X k associated with this constraint 
is strictly positive. Identifying parameter in (27) and (28) yields / {^jk,i^^j = 9 il^~+£^ , where 
/ is the function defined by (8). Replacing the value of 9k,i^f in (27) by f' 1 ( ^4^^f ) directly 
provides the following equation: 

1 - Vk,i + &i = X k gk,i F ( 9k,t ^ % — P^p] , (29) 

where F is the function defined by (15). Define Ai = {k/u kji = 0}. In other words, A\ is the set of users 
of Cell A being assigned non zero share of the band J, and Ai is the set of users of Cell A being assigned 
non zero share of IP 'a- By complementary slackness, we may write on the opposite Ai = {k/jk,i = 0} 
where E denotes the complementary set of any set E C {1, . . . K A }. After some algebra, it can be shown 
that v kji = implies fi^i = 0. Thus, 

v * e *- irkKirkA)^" 1 - <30) 

On the other hand, if v k ^ > 0, KKT conditions lead to 

vte ^ rfk F (rfk ft ) <A ' Tl <31) 

To prove that inequality (31) holds, one needs to separate the two possible cases W k ^ = and W kj i > 0. 
i) If = 0, equation (28) leads to = v k<i . Thus, (29) is equivalent to 1 — fik,i + ^i = X k g ky i, which 
implies that 1 ^ < A^ 1 since > 0. Noticing that F ( jfjjj^r Pij < 1 and multiplying this inequality by 
the previous one, we obtain the desired equation (31). ii) If W k ^ > 0, complementary slackness condition 
Hk,iW k i = along with equation (29) lead to jjL ki = = 1 + £<5j - X k g kti F ( S r^r = Br )• As function 



l-/i*,i+€<5< 

F(x) is strictly decreasing, F (if^-ft) < F (t4£tRSt) = i&- We thus obtain inec l ualit y ( 31 ) as 
well. 

To summarize, every global solution to our optimization problem can thus be characterized by the 
following set of conditions: 



1) For every k G Af. 



9k,i T-i ( 9k,i g \ _ ^-i Wk j _! ( 9k,i a \ nl > 

tt^ f \T+Wi ) ~ K ' - ^ ; (32) 



2) For every k € Ai 



9k,i 



i + W F (lTfe ft ) <A - 



(33) 



August 28, 2009 



DRAFT 



26 



3) 



\/k C k = R k , ^2"fk. 



1 = a, 



k 



1 - a 



We determine now which users are in A\ and which are in Ai- For that sake, the following conjecture 
will be revealed useful in the sequel. Define h(x) = ^^r^r ■ 

Conjecture 1. Function f(x) is strictly convex. Function h{x) is non increasing on the interval (0, 1). 

In order to validate the above conjecture, Figures 8 and 9 represent the second respectively derivative 
of / which is obviously positive, and the first derivative of h, which is obviously negative on (0, 1). We 



1.5 



Second derivative of f 



0.2 0.4 0.6 0.6 

X 



-500 



-1500 



First derivative of fi 



0.2 0.4 0.6 0.8 1 

X 



Figure 8. Second derivative of function / 



Figure 9. First derivative of function h 



show now that equations (30) and (31) are sufficient to prove that the following lemma holds. 

Lemma 1. Any global solution to Problem 1 is "binary" i.e., there exists a user L in Cell A such that 
lk,2 = 0/or closest users k < L, and 7^1 = for farthest users k > L. 

Proof: Now define L = min A2 as the closest user to the base station among all users modulating in 
the protected band Ta- By definition of L, we have 712 = • • • = 7l-i 2 = which is equivalent to the 
first part of the desired result. Now we prove the second part i.e., 7l+i,i = • • • = Jk a ,i = 0- To simplify 
notations, we define for each user k, ^,1 = By definition, L £ A2- By immediate application 
of the above KKT conditions, gi^F {9 1,2^2) = A^ 1 > gL,\F (<7z,i/3i)- As F is decreasing, we obtain 
/?2 < -^F^ 1 (gL.iPi)^. Now consider a second user k > L + l and assume by contradiction that 

k € A\. Using the same arguments, it is straightforward to show that > ^77-^ 1 (f^^ 7 (9k,iPi)j ■ 
Putting all pieces together, ^i^ 1 (f^F {h^i)) < J^ F ^ (JjlM) ■ We now prove that 

the above inequality cannot hold when k > L. To that end, we introduce the following notations. Define 
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x — 9l iPi, r = t = -r and 5 = 5 ^ . ^ • Using these notations, the above inequality reduces to 

^F" 1 (s£F (rtx)) < F- 1 {sF (x)) . (34) 

Note that in the above inequality, all variables r, s, t are strictly less than one. We now prove with the 
help of Conjecture 1 that the above inequality leads to a contradiction. In fact, Conjecture 1 states that 
function f(x) is strictly convex. As f(x) is also strictly increasing, its inverse f^ 1 is strictly concave 
strictly increasing. Therefore, for every t < 1 and for every y > 0, f~ x {ty) > tf^ 1 (y). Using the 
definition of function F(x), it is straightforward to show that the latter inequality leads to 

V(r, s, t) G (0, l) 3 , ^-F^istFitrx)) > ^-F~ l (sF{rx)) (35) 

for each real x. As function h(x) = X %-i^ is non increasing on (0, 1), it can be shown after some 
algebra [29] that function r — > ^F~ 1 (sF(rx)) is decreasing on (0, 1). As a consequence, 

V(r, s) G (0, l) 2 , ^F _1 (sF(ra;)) > F _1 (sF (x)) . (36) 

Clearly, (35) and (36) contradict inequality (34). This proves the desired lemma. ■ 
Lemma 1 establishes the "binary" property of any global solution to Problem 1. One still needs to 
prove that equations (9), (10) and (11) hold. Fortunately, these equations result directly from combining 
the above claim with equations (32) and (33). 

Appendix B 
Proof of Theorem 2 

Notations. In the sequel, x.ab represents a vector of multicell allocation parameters such that xab = 
[x A T ,x B T ] T where x A = [(W A ) T , (f A ) T ] T and x B = [(W B ) T , ( 7 B ) T ] T and where for each c = 
A,B, W c = m jl ,W^ 2 ,...,W^ 1 ,Wf CCt2 } T and 7 = [7^, 7l c ;2 , . . . , *f Ke>1 , 7 ^ 2 ] T . We respectively 
denote by Qi(x c ) = W[; l and <32(x c ) = J2k 2 tne P owers transmitted by Base Station c in the 
interference band 3 and in the protected band 5> c . When resource allocation xab is used, the total power 
transmitted by the network is equal to Q(xab) = S c Qi( x c) + <22(x c ). 

Recall that Problem 2 is nonconvex. It cannot be solved using classical convex optimization methods. 
Denote by x* AB = [x* A T ,x* B T ] T any global solution to Problem 2. 

Characterizing x* AB via single cell results. 

From x* AB we construct a new vector xab which is as well a global solution and which admits a 
"binary" form: for each Cell c, 7^ 1 = if k > L c and 7 £ 2 = if k < L c , for a certain pivot-integer L c . 
For each Cell c, vector is denned as a global solution to the single cell allocation Problem 1 when 



August 28, 2009 



DRAFT 



28 



a) the admissible nuisance constraint Q is set to Q = Qi(x^), 

b) the gain-to-interference-plus-noise-ratio in band J is set to g^i = g^^ (Qi(x^)). 

Vector xb is defined similarly, by simply exchanging A and B in the above definition. Denote by 
~x.ab = [xa T ,x# t ] t the resource allocation obtained by the above procedure. The following Lemma 
holds. 

Lemma 2. Resource allocation parameters and n* AB coincide: = x^ B . 

Proof: It is straightforward to show that jz^b is a feasible point for the joint multicell Problem 2 
in the sense that constraints C1-C5 of Problem 2 are met. This is the consequence of the low nuisance 
constraint Qi(x c ) < Qi(x*) which ensures that the interference which is produced by each base station 
when using the new allocation tlab is no bigger than the interference produced when the initial allocation 
x* AB is used. Second, it is straightforward to show that x^s is a global solution to the multicell Problem 2. 
Indeed, the power Qi(x c ) + C}2(x c ) spent by Base Station c is necessarily less than the initial power 
Qi(x*) + Q2(x*) by definition of the minimization Problem 1. Thus Q(xab) < Q{yL* AB ). Of course, as 
x. AB has been chosen itself as a global minimum of Q, the latter inequality should hold with equality: 
Q(xab) = Q( x ab)- Therefore, x.* AB and x^b are both global solutions to the multicell Problem 2. As 
an immediate consequence, inequality Qi(x c ) + C}2(x c ) < Qi(x*) + Q2{** C ) holds with equality in both 
Cells c: 

Qi(xc) + Q 2 (x c ) = Qi(x^) + Q 2 (K) • (37) 

Clearly, x^ is a feasible point for Problem 1 when setting constant Q = Qi(x^) and g^i = gt : i (Qi(x B )). 
Indeed constraint C6 is equivalent to Q\ (x* A ) < Q and is trivially met (with equality) by definition of Q. 
Since the objective function Qi(x^) + Q2(~x* A ) coincides with the global minimum as indicated by (37), 
x^ is a global minimum for the single cell Problem 1. By Theorem 1, this single cell problem admits a 
unique global minimum x^. Therefore, x^ = x^. By similar arguments, x^ = x#. ■ 
Using the above Lemma along with Theorem 1, we conclude that any global solution x* AB to the 
joint multicell Problem 2 satisfies equations (19), (20) and (21), where parameters L c , /3f, (3%, £ c for 
c = A, B in the latter equations can be defined as in Appendix III-C using values g^ i = g^ \ (Qi(x^)) 
and Q = Qi(x*). The proof of Theorem 2 is thus complete. 
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